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Abstract 

We give an exponential lower bound for the Graver complexity of the incidence matrix 
of a complete bipartite graph of arbitrary size. Our result is a generalization of the result 
by Berstein and Onn 2] for the complete bipartite graph K^ r^ f > 3. 

Keywords and phrases: algebraic statistics, contingency table, three-way transportation pro- 
gram. 

1 Introduction and the main result 

The Graver complexity of an integer matrix is currently actively investigated for its importance 
to integer programming, algebraic statistics and other applications ([3], [1], [H], [H). In particu- 
lar, from the universality of the three-way transportation program to general integer programs 
(De Loera and Onn [1]), the Graver complexity of the incidence matrix of the complete bi- 
partite graph i^3^r is particularly important. Berstein and Onn |2] proved that the Graver 
complexity g{r) for the incidence matrix of i^s^r, r > 3, is bounded below as g{r) = Vt{2'^), 
where g{r) > 17 • 2'"'^ - 7. It is a natural question to generalize this result to the complete 
bipartite graph Kt^r of arbitrary size t,r. We prove that the Graver complexity for Kt^r is 
— 1)*^), where t > 4 is fixed and r diverges to infinity. For proving our result, we employ 
double induction on r and t starting from the result of [2]. 

Let At^r denote the incidence matrix of the complete bipartite graph Kt^r and let g{At^r) de- 
note its Graver complexity. Here we state our main theorem. Relevant notations and definitions 
will be given in the next section. 

Theorem 1.1. The Graver complexity of At,r for any A<t<r is bounded from below as 

g{A,r) > {t - iy-\bt + ^ 



t-2' t 
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where 



bt = (t-2)! 15 + 



t-4 

E 

i=l 



i + 4 



We give a proof of this theorem in Section [3] after giving necessary definitions and reviewing 
relevant known results in Section [2l We conclude the paper with some discussion in Section [H 



2 Preliminaries 

In this section we summarize our notation and review relevant known results on the Graver 
complexity following Berstein and Onn Q]- 

The integer kernel of an s x t integer matrix A is denoted by kerz(^) = {x G Z* | Ax = 0}. 
Define a partial order C on Z*, which extends the coordinate- wise order < on Z^, as follows: 
For two vectors u, u G Z*, u ^ v \i \ui\ < \vi\ and uivi > for i = 1, . . . ,t. The Graver basis 
GiA) of A is the finite set of C-minimal elements in the set ker2(j4) \ {0}. 

For any fixed positive integer h, write an /it-dimensional integer vector x € Z'^* as x = 
(x^, . . . ,x^) with each block x* belonging to Z*. The type of x = (x^, . . . ,x^) is the number 
type(x) := | x* ^ 0}) of nonzero blocks of x. The h-th Lawrence lifting of an s x t matrix 

A is the following (t + hs) x /it matrix, with It denoting the t x t identity matrix: 



( ^ 







A 


h 








It 



\ 



A 
It 



(1) 



The Graver complexity of j4 is defined as 



g{A) = sup {0} U <^ type(x) 



h>l 




(2) 



Let denote the Graver basis of a matrix whose columns are the elements of Q{A) 

ordered arbitrarily. The following result shows that the Graver complexity of A is determined 
by Q{Q{A)). 

Proposition 2.1. 0/ The Graver complexity of A satisfies 

=max{||x||i : x G G{G{A))], 
where \\ ■ ||i denotes the 1-norm of a vector. 

A circuit of an integer matrix yl is a nonzero integer vector x G ker^ {A) , that has inclusion- 
minimal support with respect to kerz(vl), and whose nonzero entries are relatively prime. Let 
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C{A) denote the set of circuits of a matrix A. Then C{A) C Q{A) (cf. [7]). An integer relation 
h = {hi, . . . , hk) on integer vectors f ^, . . . , v'^ G Z* 

= hiv^ H h hkv'' 

is primitive ii hi, . . . ,hk are relatively prime positive integers and no /c — 1 of the {w*}^^;^ satisfy 
any nontrivial linear relation. By C{A) C G{A) and Proposition 12.11 we have the following 
result. 

Proposition 2.2. f^] Suppose that h is a primitive relation on some set of circuits {x^}^^i of 
an integer matrix A. Then the Graver complexity of A satisfies g{A) > Yli=i ^i- 

In this paper we consider the Graver complexity of the incidence matrix At^r for the complete 
bipartite graph Kt^r- Let 1^ = (1, 1, . . . , 1) denote the 1 x t matrix consisting of I's. Then the r- 

(r) 

th Lawrence lifting At^r = 1* of It is the incidence matrix of Kt^r- In algebraic statistics, At^r is 
the design matrix specifying the row sums and the column sums of a two-way contingency table. 
Another Lawrence lifting (1^^^)^'^) of 1^^^ is the design matrix for no-three-factor interaction 
model for t x r x h three-way contingency tables It is also the coefficient matrix 

for the three-way transportation program. The Graver complexity g{At^r) = di'^t^^) gives the 
bound of complexity of the Graver basis for the toric ideal associated with the no-three-factor 
interaction model for t x r x h three-way contingency tables as /i — >■ oo. 

We employ below the following notation, where t, r are positive integers. Let 

V := {vi, . . . ,vt}, U := {ui, . . . ,Ur}. (3) 

Then V (BU and V xU denote the set of vertices and the set of edges of the complete bipartite 
graph Kt^r, respectively. They index the rows and the columns of the incidence matrix Ai ^ 
of Kt^r- Here we explain interpretations of a circuit of At^r referring to We interpret each 
vector X E Z^^^ as: 

1. an integer valued function on the set of edges V x U; 

2. a t X r matrix with its rows and columns indexed by V and U. 
With these interpretations, x is in C{At^r) if and only if: 

1. as a function onVxU with the following properties: its support is a circuit of Kt^r, along 
which its values ±1 alternate. It can be expressed by the sequence (wj^ , Ui-^ ,Vi2,Ui^, . . . , , Uj, ) 
of vertices of the circuit of Kt^r on which it is supported, with the convention that its 
value is +1 on the first edge {vi-^,Ui-^). 

2. as a nonzero matrix with the following properties: its elements are 0, ±1, its row sums 
and columns sums are zeros, and it has an inclusion-minimal support with respect to 
these properties. 

The following example is the base case for our inductive argument for the lower bound of 
the Graver complexity. 
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Table 1: Circuits for ^3^4 in a 3 x 4 matrix form 
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V2 











-1 


1 
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x^ 


= (fi,n4,f2 


U2,V3,Ui) = 
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-1 


V2 
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-1 








V3 
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-1 
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x^ 


= (fl,Ul,f2 


U2,V3,U3) = 


-1 


1 








V2 











-1 
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-1 
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Vi 


x^ 


= (fl,U3,f2 


U4,V3,U2) = 








-1 
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V2 











1 





-1 
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1 





-1 
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= {vi,U2,V2 







-1 


1 





V2 














-1 


1 


V3 



Example 2.1. J^/ Let t = 3 and r = 4. Consider seven circuits {x*}-'^]^ in Table[l\ (written 
in a 3 X 4 matrix form) of ^3^4 = (1, 1, 1)*^^^ . They satisfy a primitive relation 

x^ + 2x^ + 3x^ + 3x^ + 5x^ + + 7x'^ = 0. 

Therefore from Proposition \2.S\ 

5(^3,4) >l + 2 + 3 + 3 + 5 + 6 + 7 = 27. 

3 Proof of the main theorem 

In this section we give a proof of our main theorem. Our proof is based on recursive construction 
of primitive relations for circuits of At^r- We need recursions for t and r, separately. In Lemma 
13.11 we give a recursion for t and in Lemma 13.21 we give a recursion for r. 

Lemma 3.1. Let t > 4. Suppose that there are circuits {x*}^^;^ of At^t+i = ll*^^"* admitting a 
primitive relation h, where the k-th circuit and the k-th coefficient are 

X^ = {vi,U2,V2,U3, . . .,Vt,Ut+l), 
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hk = 1. 



Then there are circuits {a;*}^^^ of At+i^t+i = '^t+i admitting a primitive relation h, where the 
{k + t)-th circuit and the {k + t)-th coefficient are 

x'^+'t = {vi,Ul,V2,U2, ■ ■ ■ ,Vt+l,Ut+l), 

hk+t = 1- 

Proof. Using the natural embedding of Kt^t+i into i^Tf+i^t+i, we can interpret circuits of the 
former also as circuits of the latter. Put 

y' = x\ Vi = i,...,fc-i, 

and define 

yfc+j-i ^ [vi,Uj,vt+i,Uj+i), Vj = 1, . . . ,t, 

y''+^ = {vi,U2,V2,U3, . . . , Vt, Ut+1, Vt+1 , Ui) . 

Table [2] displays {y'^+^-H*! \ as matrices, where 

p= (1,0,..., 0,-1)^ GZ*+^ 
Blank entries are zeros. Note that these circuits satisfy X^jii 2/*^^"' ""^ ~ 



Table 2: Circuits for recursion on t 





Ui 


U2 


U3 


Ut Ut+l 


= {vi,Ui,Vt+l,U2) = 


P 


-p 






y/^+l = {vi,U2,Vt+l,U3) = 




p 


-p 








yk+t-1 ^ {vi,Ut,Vt+l,Ut+l) = 


p —p 


yk+t-1 ^ [vi,Ut,Vt+l,Ut+l) = 


p —p 




-1 


1 
-1 


1 




yk+t ^ (^v^^u2,V2,U3, . . . ,Vt,Ut+l,Vt+l,Ui) = 


1 






-1 1 

-1 1 

-1 



Suppose that h € Z'^^* satisfies 

hi = hi, Vz = 1, . . . , /c - 1, 
hk+j-i = hk, yj = 1, . . . ,t + 1. 
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Then 

k+t k-1 t+1 k-1 

i=l i=l j=l 1=1 

Therefore h is an integer relation of circuits {y'jfif . 

Next, we show that h is primitive. Suppose that h' G is a nontrivial relation on 

the {y^}i^i- Without loss of generality we may assume that the are relatively prime 

integers, at least one of which is positive. We look at the row of vt+i- Then it follows that 

h'k = h'k+l = • • • = /ife+f 

Therefore 

k+t k-1 t+1 k-1 

= J2 Ky' = E + E y'^''' = E ^^^^ + ^'k^'- 
1=1 1=1 j=i 1=1 

This is an integer relation on and because h is primitive, 

h^ = hi, V^ = 1, . . . , fc. 

Therefore h' = h and /i is primitive. 

Now apply to a permutation of columns so that y'^^* becomes {vi,ui,V2,U2, ■ ■ ■ , vt+i,ut+i). 

For i = 1,. .. ,k + t, let he the circuit of At+i^^+i which is the image of under this permuta- 
tion. Then also satisfy the primitive relation Yli=i ^i^' = with the same coefficients 
h. This completes the proof. □ 

Lemma 3.2. Let r >t> 4. Suppose that there are circuits {x^}^^^ of Af^r = ^i^^ admitting a 
primitive relation h, where the k-th circuit and the k-th coefficient are 

X^ = {vi,Ur-t+l,V2,Ur-t+2, ■ ■ ■,Vt,Ur), 

hk = 1. 

Then there are circuits {x^}^^l~^ of At^r+i = 1*'^'''^'* admitting primitive relation h, where the 
{k + t — l)-th circuit is 

^fc+t-l ^ {vi,Ur-t+2,V2,Ur-t+3, ■ ■ ■ ,Vt,Ur+l) 

and the elements of h are 

hi = {t- l)hi, yi = l,...,k-l, hk = hk+i = ■■■ = hk+t-i = hk = I. 

Proof. Using the natural embedding of Kt^r into Kf^r+i, we can interpret circuits of the former 
also as circuits of the latter. Put 

yi = x\ Vz = l,...,fc-1, 
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Table 3: Circuits for recursion on r 





■ ■ ■ Ur^t-\ 


-1 U-r-t^ 


-2 «r-t+3 • 




y'' = {vi,Ur+i,V2,Ur-t+2, - ■ ■ 


,Vt-l,Ur-l,Vt,Ur) = 









. g* g^ 


y*=+l = [vi , Ur-t+l , V2 , U^+i , . . 


. ,Vt-l,Ur-l,Vt,Ur) = 


... 







. q^ g^ 








■ ■ ■ ,Vt-l,Ur+l,Vt,Ur) = 


... 


l' 


q"^ 


. q^ q^ ^ 


yk+t-l ^ {vi,Ur-t+l,V2,Ur-t-^2,- 




... 




q'^ 


. g* 



and for all j = 1, . . . , t, let y^+^~^ denote vectors obtained by changing vertex Ur-j+i of to 
Ur+i. Table El displays these circuits as matrices. Here for each i = 1, . . . , i, g* G Z* denotes a 
vector satisfying 

i = 1, 4+1 = -1, 

and the rest are zeros. Here we identify t + 1 with 1. 
Notice that 

Define 

hi = {t — l)hi, \/i = 1, . . . ,k — 1, 
hk+j-i = hk = 1, Vj = 1, . . . , t + 1. 

Then 

k+t-l fe-1 t k-l 

hy' = Y^it - i)hiy' + Y hky'^'-' = Y + hkx' = 0. 

1=1 i=l j=l 1=1 

Therefore h is an integer relation on circuits {y^}i^i~^- 

Next we show that h is primitive. Suppose that h' G Z'^"'"*"^ is a nontrivial relation on 
the {y^}i^i~^- Without loss of generality we may assume that {h[}^^^~^ are relatively prime 
integers, at least one of which is positive. Consider the column of Ur+i- Then 

^fe = ^fc+i = ■ ■ • = 

Therefore 

k+t-l k-1 t k-1 

= E ^'^y' = E h'^y' + ^'^ E y'^'~' = E + - 1)^'^^'- 

i=l 1=1 j=l 1=1 

This is an integer relation on {x'}^^^. Therefore there exists a G Z such that 

h'i = ahi, Vi = 1, . . . , /c — 1, (4) 
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{t — = ahk = a. 



(5) 



Since hi > for all i and there is an i such that h'- > 0, equations (jH and ([5]) imply a > 0. 
Therefore ^ and ^ imply that h'- > for all i. Hence h is primitive. 

Now apply to a permutation of columns so that y^+*~i becomes (f i, nr_t+2) ^'2) ^^r-t+s, 

. . . , ft, Ur+i). For i = 1, . . . , A; + t — 1, let X* be the circuit of At^r+i which is the image of 
under this permutation. Then {x^}^^f~^ also satisfy the primitive relation X]f=i~^ ^j^* ~ ^ 
with the same coefficients h. This completes the proof. □ 

We are now ready to prove Theorem 11.11 In the proof we use the following notation. Let 
■^{{x^}i=i) = {^*}f=i ^'iid ^{h) = h = {hi, . . . ,hk+t~iA) denote circuits of At+i^t+i and 
the primitive relation which are obtained by the operation of Lemma l3.ll to circuits {x*}^^]^ 
of Af^t+i and the primitive relation h. Note that = \\h\\i + t. Furthermore let 

=c/'({x*}JLi) = {x^lfi"*"^ and ^'{h) = h = (hi,..., hk+t-2, 1) denote circuits of At^r+i and the 
primitive relation which are obtained by the operation of Lemma [3.21 to circuits {x*}^^^ of At^r 
and the primitive relation h. Note that ||^'(/i)||i = {t - l){\\h\\i -l)+t. 

Our proof uses induction on t,r. We will construct a primitive relation /i^*^*^) on circuits 
^{txr) of ■^t,r by induction. Therefore we obtain g{At^r) > H''-^*^'"^ ||i- Our induction is illus- 
trated in Figure [TJ There, a down arrow corresponds to the operation of Lemma 13. H and a 
right arrow corresponds to the operation of Lemma |3.2[ 

||/l(3x4)||^ ^ ||/l(3x5)||i ^ ||/i(3x6)||i ^ ... 
i 

i 

||/l(5x5)||^ ^ ^ ••• 



Figure 1: Induction on t,r 



Proof of Theorem \l.l\ By induction on t we will prove that for all t > 4 there exist k{t) 

^{txt) = {^\tyct) 
kit) 



— 2t + 2 circuits ^uxt) = {^)txt)^'i=i ^ C{^t,t) and the primitive relation /i(*x<) such that 



^(ixt) = {vi,Ui,V2,U2,. ■ ■ ,Vt,Ut), 

kit) 



1=1 



, (txt) _ . 

'''kit) ~ 
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Exchange and of circuits of Example 12.11 and apply to the circuits a permutation of 
vertices so that 

= {vi,U2,V2,U3,V3,U4). 

Let =^(3x4) = {x'^(^^^^j}J^i be the image of {x*}J^]^ under this permutation. The primitive 
relation h^^^^^ on these circuits satisfy 



^(3x4) ^ (7,2,3,3,5,6,1). 



Notice that hi^^^^ = 1 holds. 



Let Jr(4x4) = ■«^(^(3x4)) and /i^^^^) = ^(/i^^x^)) denote the image of Jr(3x4) and /i^^x^) 
under the operation of Lemma 13.11 Then we have = (7, 2, 3, 3, 5, 6, 1, 1, 1, 1) E Z^^ and 



2^(4x4) = {vi,Ui,V2,U2,V3,U3,V4,U4), 

'ho — 

|[/^(4X4)||^ = ||/l(3x4)||^ + 3 = 30. 

Therefore we have verified the initial condition at t = 4 for the induction. 

Suppose now that the result holds for i > 4. Let ^(tx{t+i)) = -^'{•^(txt)) and /i(*x(i+i)) _ 
denote the image of ^[txt) and /i(*xt) ^ iJ^i^) under the operation of Lemma 13.21 

k{t)+t-l , N 

^(t X (t+1)) = (^'l , ^^2 , 2 , Us , . . . , , ) , 
, (tx(t+l)) _ . 

||/l(*>^(*+l))||l = (t-l)(||/l(*>^*'||l-l)+t 

follows from LemmaO Now let ^((j+i)x(m)) = ^(^(tx{t+i))) and = ^(/i(*x(*+i))) 

denote the image of <^(^x(t+i)) and /i(*x(t+i)) under the operation of Lemma |3. 11 Then 

^((t+l)x{t+l)) = (^^l,Ul,V2,U2, . . . ,?;t+l,ttt+i), 

, ({t+l)x(t+l)) _ 
"'fc(t)+2t-l ~ ^' 

||^((m)x{m)) = (t _ i)(||/i{*xt) _ 1) + 2t 

/ (m)-4 
= ((i + l)-2)! 15+ 



follows from Lemma O Here k{t + 1) = + 2t - 1 and A:(4) = 10 imply kit) = _ 2t + 2. 
Therefore the result holds for t + 1. Henceforth, let ht = 
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We fix t > 4 arbitrarily. We prove by induction on r that, for all r > t, there are circuits 

(txr) = {^ltxr)^i 



'^(txr) = {^}txr)'^i=i *- ^(.^t,r) and the primitive relation /i^*^*") such that 



■^(txr) 


= {vi,Ur-t+l,V2,U, 


m 

2_^"'i ^{txr) 
i=l 


= 0, 




= 1, 




= (t - ly-' (h + 



t - 2 



The initial condition of the induction, at r = t, follows from = bt- 

Suppose now that the result holds for some r > t. Let ^(^tx{r+i)) = -^'{^(txr)) and 
^(tx(r+i)) _ (jg' denote the image of ^(^txr) and h^^^^^ under the operation of Lemma 
[321 Then 

k(t) + t-l , X 

^(ix(r+l)) = \y\,Ur-t^2,V2,Ur-t+Z, ■ ■ ■ ,Vt,UrJ^\), 

, (tx(r+l)) _ . 
"■k{t)+t-l ~ 
||/j(tx(r+l)) = (t _ l)(||/i(*xr) _ 1) + t 

= (^-i)((t-ir-(5. + ^)-^-i)+* 
= - ir f 6, + ^ ^ ^ 



t - 2y t - 2 

follows from Lemma 13. 2i Therefore the result holds for r + 1 and 

r-tr 1 ^ 1 



g{At,r)>{t-lY-\bt + 



t-2' t-2 

follows from Lemma |2.2[ □ 



4 Discussion 

In this paper we provided a lower bound in Theorem 11.11 by the induction on t, r. Here we 
discuss some ideas for improving our lower bound. 

Look at Figure [J again. On the step H/i*-^^^-* ||i — > ||i, we can construct a larger 

primitive relation than the relation constructed in the proof. 

Example 4.1. Let {x'^Yi=i denote the circuits in Example \2.1[ Using the natural embedding 
of K3^4 into Ki^i, let 

= x\ V« = 1,...,6. 

and for i = 7, . . . , 10, we define as shown in Table^ Then {x'^}]^-^ are circuits of hdiii'^^f^) 
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Table 4: Circuits of 
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U2 




M4 















1 





-1 


^ 1 




77/1 7/1 


7 Jo 


7/ /I 1 — 














Vo 










-1 
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V-\ 

^0 










1 


-1 








^4 






















V^ 
^1 




n3,f4, 









-1 


1 





V2 






















V3 













1 


-1 





va 






















Vl 


X? = (t>3, 


n4,f4, 


















V2 
















-1 


1 


















1 


-1 


va 













1 





-1 


Vl 


= (fi,n2,f2 






,^4,^4) = 





-1 


1 





V2 










1 





-1 





V3 










-1 








1 


VA 



and h = (2, 4, 6, 6, 10, 12, 7, 7, 7, 7) is its primitive relation. Then, by Proposition \2.^ 

g{lt^) >2 + 4 + 6 + 6 + 10 + 12 + 7 + 7 + 7 + 7 = 68. (6) 
Equation ^ is sharper than the evaluation in Theorem \l.l\ 

We could start induction from circuits and its primitive relation in Example I4.1i Then we 
obtain a sharper evaluation for some small t, r. However, unfortunately it turns out that, if 
we start from Example 14.11 then on the step ||i — > Ij/i'-^^^''^"'^^^ ||i, we can not obtain an 

exponential lower bound. Therefore we did not use Example 14.11 in the proof of Theorem 11.11 
However this example suggests that there may be some other better initial set of circuits for 
our induction. 
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